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ABSTRACT 


Several radiation models for low density nonequilibrium 
hypersonic flow are studied in this report. It is proposed that 
these models should be tested by the 3-D VRFL code developed at 
NASA/JSC. A modified and optimized radiation, model may be 
obtained from the testing. Then, the current VRFL code could be 
expanded to solve hypersonic flow problems with nonequilibrium 
therma 1 radiation. 

INTRODUCTION 


When hypersonic flight vehicles, such as Aeroassisted Qrbiter 
Transfer Vehicles (AOTVS), traveling in low density env 1 r o nmen t , 
strong bow shocks occur adjacent to the vehicles. Behind the 
shock front, the flow involved could be at thermal and chemical 
nonequilibrium. At the extreme high temperatures, the effect of 
nonequilibrium gas radiation may be prominent. Determining the 
magnitude of radiation in this environment is considered to be one 
of the most important problems associated with AOTVS (1-3). 

There are many research works dealing with thermaL and 
chemical honequ i 1 i br ium hypersonic flow, some of the recent papers 
are given in the references (4-10). Although the nonequilibrium 
radiation is part of the nonequilibrium precesses, it is usually 
assumed as equilibrium radiation or neglected completely in the 
analysis. The relatively few work in nonequilibrium radiation 
research will be mentioned in the later sections. 

The purpose of this report is to give a brief discussion of 
nonequilibrium radiation in low density hypersonic flow, with 
various models for radiative interaction, numerical solution 
methods, and some of the available codes. In the next section, 
the solution methods used for the general nonequilibrium 
hypersonic flow will be discussed, this is followed by the section 
dealing with radiation. 

THERMAL AND CHEMICAL NONEQU I L I BR I UM 
HYPERSONIC FLOW 

There are several numerical methods which can be used for the 
nonequilibrium hypersonic flow (11). The Direct Simulation Monte 
Carlo method (12,13), based on probab i 1 ist ic collision model, 
requires excessive amount of computer time to solve a two - 
dimensional flow problem so it may be not able to solve a three - 
dimensional flow problem. The Viscous shock Layer method (14, 
15), which used the simplified conservation equations in the shock 
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layer, fails to predict the flow for the whole vehicle* To solve 
the complete flow field around a hypersonic vehicle, the full 
Navier stokes equations should be used,, 

Full Navier-Stokes Equations 

In order to calculate radiation from the nonequi 1 ibr ium 
region of the flow, one needs to know the thermal and chemical 
state of th gas in the region- This includes the nonequ i 1 i br i urn 
species concentrations and various temperatures to express the 
different energy modes of the species- These state variables can 
be calculated by simultaneously solving a set of mass, momentum, 
and energy conservation equations, the so— called full Navier - 
Stoke equat ions (2,16,17). The equations consist of one mass 
conservation equation for each species, one overall momentum 
equation, and several energy equations. The number of energy 
equations should equal to the number of temperatures which 
characterize the various energy modes of the species. For 
example, when using a two-temper a ture model of one electron 
temperature and one heavy-par t i c I es temperature, the two energy 
equations could be one electron energy conservation equation and 
one overall energy conservation equation. If the two- temper a ture 
model used involving one electron-vibrational tempera ture and one 
rotational-translational temperature, the two energy equations 
could be one e 1 ec ton-v i brat i ona 1 energy conservation equation and 
one overall energy conservation equation. For a three- tempera ture 
model of an electron temperature, a vibrational temperature for 
all species, and a t r ansi a t i ona 1 -r o t at i ona 1 temperature for all 
species, then there are three energy conservation equations for 
them. In the literature, the thermal nonequ i 1 i br i urn model could 
involve more than three temperatures (18), then the computational 
time has to increase greatly. 

In the overall energy equation, there is a term represents 
the contribution of thermal radiation. This term is expressed by 
a divergence of a radiation heat flux vector. To evaluate this 
radiation heat flux vector, a radiative transport equation is 
needed. The coupling of these equations makes the nonequ i 1 i br i urn 
flow problem very difficult to solve. The radiative transport 
equation is discussed in the later section. 

Rate Equations 

To formulate the full Navier-Stokes equations for 
nonequilibrium flow, various rate expressions for source functions 
should be specified. In the mass conservation equation for each 
species, the source term is the mass production rate of species 
due to chemical reactions. The chemical reaction rate is a 
function of the various temperatures used for species. In the 
energy equation involving translational - vibrational energy 
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exchange, the relaxation time in the rate equation is of the 
Landau - Teller type <2, 17). For* the el.ectron-vibrat ional energy 
exchange, the relaxation time expression is a modification of the 
Landau-Teller expression (2,19). 

NONEQUILIBRIUM GAS RADIATION 


As mentioned in the earlier section, to include a radiation 
energy term in the overall energy equation, an additional equation 
is needed, this is the equation of radiative transfer. The 
coupling between this radiative transfer equation and the energy 
equation in gas dynamics is of a complicated i n tegr a 1 -d i f f er en t i a l 
form. When nonequilibrium requirement is included, it adds more 
complexity to the problem. In the following sections, the basic 
equation, various simplication and solution methods are briefly 
d i scussed . 

The Equation of Radiative Transfer 

The equation of radiative transfer is a continuity equation 
for the number density of photons in a specific solid angle and at 
a specific frequency (17). The equation can be formulated in 
terms of intensity, emission and absorption coefficients. In a 
nonequilibrium environment, a knowledge of the emission and 
absorption coefficients is a requisite to a careful solution of 
the equation of transfer. The discussion of various spectral 
models for emission and absorption coefficients is in the next 
section. since the radiation transport is coupled with flow 
equations, for a formal solution procedure, one must rely on an 
iterative method to solve the system of flow equations with 
radiation tr anspor t . 

Radiation Transport Models 

There are a number of spectral models for radiative 
interaction computations in the literature. These models are 
briefly outlined in the following: 

Line By Line Transport Calculation Park, etc (3,20,21) has 
developed a model to calculate radiative properties for 
nonequilibrium air. the model assumes that the radiation 
transport is not strongly coupled with the flow equations, so the 
heavy-particle translational temperature, vibrational teniperature, 
electron temperature, and species mole fractions are solved 
independently from flow equations. These items will be used as 
inputs to calculate the number densities of internal state of 
atoms and molecules by using the so called quasi steady-state 
assumption. A detailed line-by-line calculation is performed to 
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evaluate emission and absorption characteristics independently, 
using the calculated, nonequilibrium, excited internal state 
density values. This model requires a tremendous amount of 
computer time and memory space for the required computations. It 
becomes even more so when four vibrational temperatures were used 
for the gas (18). It is very difficult to apply this model to a 
general 3-D hypersonic flow problem. 

The Step-Function Absorption Coefficient Model A step function 
model has been developed to compute the radiation absorption 
coefficients by Zoby, Sutton, etc. (22). In general, the model 
includes, free-free, bound-free, and bound-bound transitions of 
various atoms and molecules. The inputs to the model are species 
niimber density and temperature of the flow field. The frequency 
dependence of the absorption coefficient is represented by a set 
of steps, these steps has a fixed, but not necessarily equal 
widths chosen to resolve the detailed spectra adequately. The 
total absorption coefficient of a particular step is a summation 
of all the average absorption coefficient of individual energy 
transitions in that step. The average absorption coefficient of 
an individual transition can be obtained by integration the 
absorption coefficient over frequency in that step. When the step 
width is properly chosen, there should be analytic expressions to 
approximate the absorption coefficient in that step for the 
purpose of integration. Thus the frequency dependence for each 
step- becomes a set value, and the total absorption coefficient for 
that step can now be expressed in terms of species number density 
and temper a ture . When this is done, the radiative transfer 
equation is well expressed in each frequency step. In Reference 
EE, a 58 steps were used to cover the frequency range from 0 to 
1 7 ev . 


Although this model is mainly applied to equilibrium 
radiation, it includes the frequency dependent absorption 
coefficient (22,23,24). In addition, with some mod i f i ca t i ons , 
this method can be applied to nonequilibrium radiation (25). 

There is a similiar but simpler approximation called the band 
approximation model for absorption coefficient (26,27). This 
model is also for equilibrium raidation and usually is a two-band 
approximation model . 

Parametric R a diation Model The nonequilibrium radiation 
environment corresponds to conditions where population densities 
of the energy levels deviate from the equal ibrium Boltzmann 
distributions based on a single temperature. At sufficiently low 
densities, the probability of a radiative transition becomes 
comparable with the probability of a correspond i ng collisional 
transition. Unless the gas is optically thick, the emission of a 
photon is not balanced by its inverse. Consequently, the 
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population distribution among the energy levels departs from the 
predicted by the Boltzmann equation (25,28). In the formulation 
of the radiative transfer equation, Tiwari, etc. (25,29) devised a 
model that the absorption coefficient of different species is 
calculated by using the 58-step absorption model developed by 
Sutton (22), while the major nonequ i 1 i br i urn effect will enter 
through the source function in the equation. For a multilevel 
energy transitions model, this is done by introducing a parameter 
n which represents the ratio of the collisional relaxation time 
and the radiative lifetime of the first excited state, and the 
higher level energy transitions can be related to this parameter. 
The nonequilibrium radiation become important for conditions when 
n = 0 (1). The key to this model is the determination of 
collisional relaxation time between various particles and 
radiative lifetime of the excited states. Additional comparison 
and analysis are needed for this model. 

Other Radiation Models There are a number of other radiation 
models in the literatures (30-32). Further investigations and 
testings are needed for these models. In the field of 
nonequilibrium hypersonic flow with thermal radiation, it is 
desirable to have more • theoret ical analysis, computational code 
developments, and experimental data collections. 

Integration of Transfer Equation 

.To obtain the radiant heat flux vector in the energy 
equation, the radiative transfer equation needed to be integrated 
spectrally as well as spatially. This is very difficult to 
perform on a flow field over a whole vehicle. For an 
approximation, the tangent slab approximation technique is 
commonly used (20,22,23,29). This approximation implies that the 
radiative energy transfer along the body is negligible in 
comparison to that transferred in the direction normal to the 
body. Other approximation such as optical thin or/and optical 
thick assumptions may also be used to solve the problem (18). In 
general, additional works are needed in the area of fast spatial 
integration of transfer equation. 

CONCLUSIONS 


Several radiation models were found in the literature servey. 

These models need to be tested by some existing 3-D Nov i er -s to kes 
computer codes. It is proposed that the VRFL code is used for 
this purpose (7). From the result of testing, comparisons can be 
made on accuracy of prediction, requirement of computer time and 
memory space, and generality of the model. To save computer time 
and to facilitate easy implementation, these models should be 
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simplified but retained their essential features. In the computer 
tests, the flow field may be limited to the stagnation region of a 
symmetric blunt body for tangent slab approx imat ion, use minimum 
number of temperatures for thermal model, and assume radiation 
heat transfer equation and gas dynamic equations are uncoupled. 
From this testing and additional analysis, a modified and 
optimized radiation model may be obtained. 

After the optimized radiation model is incorporated into the 
VRFL code, the heat transfer rate to a entire hypersonic vehicle 
should be calculated. The radiation heat transfer equation and 
gas dynamic equations are still uncoupled but the tangent slab 
approx imat ion is removed. The purpose of this test is to find a 
numerical model for the spatial integration of the radiative 
tr ansfer equation. 

Finally, the expanded VRFL code can he applied to vehicles 
such as AOTVs. For these computations, the coupling of radiation 
and gas dynamics should be maintained. The purpose of calculation 
is not only to accurately predict the heat transfer to the 
vehicle, tout also to find the effect of radiation to the 
hypersonic flow field. 

New radiation model and numerical solution technique may 
still need to be developed. At present there are still too many 
unanswered questions in the field of nonequilibrium hypersonic 
flow with thermal radiation. 
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